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ABSTRACT 


Recently the dispersion curve of the elementary excitations in Helium II has been studied 
by observing the angular and energy distribution of cold neutrons scattered from liquid 
helium, in Stockholm [1,2] in the momentum region 1.50<%<2.14 A* at a helium tem- 
perature of 1.44°K, and at Los Alamos [3] in the region 0.55<%x< 2.36 A”* at a helium 
temperature of about 1.1°K. The latter workers also studied the variation with helium tem- 
perature up to 1.8°K of parameters of the dispersion curve such as the width of the roton 
line and the minimum excitation energy A in the roton region, where the dispersion curve 
according to Landau is given by 
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A is the minimum energy to create a roton, p the momentum of the roton and y its effective 
mass. Independently of the Los Alamos work the variation with temperature of the dispersion 


eurve has been studied at the Stockholm reactor in the temperature region 1.44< 7’<2.37°K. 
The observed temperature dependence of A is very well fitted by the functions 
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A 8.5 (fiyyeyare; T) for 1.4<7<1.9°K 
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Tweet rar r for 1.8<T<2.15°K 


The experimentally observed width 6H of the roton line is well described by 


8.5 
6B | s02-Vre? for T'<2°K 
eas pe —85)8/2 
Broth: r) for T>2°K 


Introduction 


Shortly after the theoretical prediction by Feynman and Cohen [4] that cold 
utrons could be used to excite single excitations in liquid helium, experiments 
1, 2,3] were performed which confirmed this idea. Moreover, these experiments 
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successfully showed that the basic concepts of Landau [5], Feynman [6], and 
Brueckner and Sawada [7] were correct: the energy vs. momentum relation for 
the excitations showing a marked minimum at a momentum of about 1.90 A“?, 


ix A ; 
the minimum energy ],,’ necessary to create a roton being about 8.1-8.6°K at 
B 


helium temperatures between 1.5 and 1.1°K. Theoretically wery little is known 
about what happens at higher helium temperatures, up to the A-point, when the 
density of the excitations increases fast. According to the Landau-Khalatnikov 
theory the roton ye REET N, (number of rotons per unit volume) should 


increase as V7'-e 7 with helium temperature. This would result in a decreasing 

mean life time of the rotons and consequently decreasing mean free path. This 

in turn would correspond to a roton line width increasing with temperature as 
A 


VT +e 7. This law would then be valid up to helium temperatures of about 1.8°K. 
In the neighbourhood of the A-point there is no theoretical prediction. Also there 
is no theory for any variation of the “gap” energy A with helium temperature, 
except the attempts of Brueckner and Sawada to calculate the dispersion curve 
for a dense Bose gas. For these reasons we undertook a systematic study of the 
temperature variation of the dispersion curve using the cold neutron technique. 


Experiments 


a) Experimental method and resolution of neutron spectrometer 


In earlier papers [1,2] the principles used to produce single excitations in 
Helium II by cold neutrons have been described in detail. In the present inves- 
tigation the same apparatus [8] and the same technique of measurement have 
been used. The angular and energy distribution of the beryllium filtered neutron 
beam scattered from a liquid helium sample is analyzed by means of a slow 
chopper time-of-flight spectrometer. The neutron spectrum impinging upon the 
helium sample is characterized by a sharp edge at a neutron wave length of 
3.95 A corresponding to the last Bragg cut-off of beryllium. The resolution width 
of the spectrometer is such that the full width, (6¢)), at half maximum of the 
resolution function, assumed gaussian in shape, is 100 +5usec, when the time-of- 
flight of 4A neutrons is about 3000usec. The primary spectrum scattered elas- 
tically from a vanadium sample is given in Fig. 1, where the resolution width 
is shown as well as the position of the beryllium break. The beryllium cut-off 
very nearly corresponds to the position of the point of inflection of the sharp 
edge, which is rounded off by the finite resolution function. The determination 
of the true break position is slightly complicated by the presence of an alumi- 
nium “break” at 4.05 A. Another aluminium ‘“‘break’’ shows up at a wave length 
of 4.65 A (compare Fig. 1 and reference 8). This complication arises from con- 
structional aluminium in the reactor tank and experimental channel used. Cal- 
culation shows, however, that. the correction due to this cause is small, amounting 
to about 3ysec in 3000usec. Another complication arises from the fact that 
the filtering action of the beryllium filter of 20 cm thickness is incomplete, 
resulting in a minor transmission at the 3.58 A beryllium ‘“‘break”’ corresponding 
to the peak at a flight time 2850 usec in Fig. 1. This fact also causes a minor 


50 


ARKIV FOR FYSIK. Bd 15 nr 4 


t, «3069 psec 


Counts 
Uchennet 


i] 
' 
' 
| 
a 

($t}s100usec 


2000 2500 3000 3500 4000 4500 5000 
psec 


Fig. 1. Spectrum of beryllium filtered neutrons as observed. by scattering from vanadium. The 
half width value of the resolution function is given as the projection of the tangent at the point 
of inflection on the time-of-flight axis and is found to be 100 usec. 


uncertainty of the order 2 or 3ysec in the definition of the 3.95 A cut-off. The 
corrected value for the time-of-flight (f)) of the 3.95 A neutrons used in the 
treatment of the present data is 3069 +5ysec. The accurate calculation of the 
shape of a triangular neutron spectrum as seen after the application of a gaussian 
resolution function and a discussion of the corrections is given in the appendix. 
If the flight time of the neutrons after the scattering process in the helium 
sample is denoted by f, the laws of energy and momentum concervation give for 
the excitation energy (7K) and the corresponding momentum value (x A’): 


Ma toth y, 
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where m is the neutron mass, kg is the Boltzmann constant, A, is the beryllium 
cut-off (3.952 A at liquid nitrogen temperature), 9 is the angle of observation, 
and At=t,—t, is proportional to the wave length shift caused by the scattering 
process (compare reference 2). 

If any change in the slope of the scattered beryllium cut-off is observed, such 
that the break appears symmetric but of a width (dt),, the true roton line width 
(dt),, expressed in time units, may be calculated from 
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(6t), = V(6t)?— (66/8 (3) 


This formula is strictly valid if the resolution function as well as the roton line 
are both of gaussian shape, as proved in the appendix. The roton line width, 
AT, in °K is then simply calculated from 


2 
éT =—_, 3: (6), 7 (4) 
if (6t),<t;, which is always true in the present case. 


b) Measurements 


First the arm of the spectrometer was fixed in such a position that the angle, 
6, of observation of the neutrons scattered from the helium cryostat with respect 
to the incoming neutron beam was 76.3°. From the earlier measurements [2] of 
the dispersion curve it was known that this angle of observation corresponds to 
the minimum point (A, x») of this curve at a helium temperature of 1.44°K. 
Then a series of runs were taken at helium temperatures of 1.50, 1.65, 1.73, 1.85, 
1.94,..1:98, 2:05: 2446917; and. 237. keeping the angle of observation fixed. 
A run at 4.2°K had been made earlier [2]. The experimental procedure involved 
the usual careful background measurements between every helium run, and check 
of the calibration of the time scale between the helium runs using the elastic 
scattering from vanadium. As demonstrated in Fig. 2 two effects were observed 
as the helium temperature was increased: 

1. The slope of the scattered beryllium break, characterized by (6t);, increased, 
corresponding to an increasing roton line width. 

2. The position of the break, t;, shifted continuously towards shorter flight times, 
corresponding to an apparently decreasing value of the excitation energy. 

The question then arises whether the excitation energies which can be calculated 
from these series of observations really correspond to the minimum, A, of the 
dispersion curve. To investigate the shape of the dispersion curve in the neigh- 
bourhood of the minimum a series of observations at different angles was there- 
fore made at a helium temperature of 2.03°K. The results of these measurements 
are given in Table 1. The values of the excitation energies (7K) and the cor- 
responding momenta (x A’) have been calculated using equations (1) and (2), 
taking t) = 3069 +5 usec. 


Table 1. Dispersion curve determined at a helium temperature of 2.03°K. 


Time-of-flight 


Angle of Time shift of Excitation Momentum 
‘ for scattered fe 
observation break position energy value 
neutrons 3 = 
At psec PENS xA 
ty sec 

90.0 3320+ 10 251412 8.84 0.4 2.164 0.5 
83.0 3266 + 10 197+ 12 7140.4 2.044 0.5 
80.0 3255+ 10 186+ 12 6.74 0.4 1.98 + 0.5 
76.3 3244+ 10 175+ 12 6.4+ 0.4 1.91+0.5 
74.0 3263+ 15 194+ 16 7.0+0.5 1.86 +0.5 
71.5 3278 +15 209 + 16 7.5+0.5 1.80+ 0.5 
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Fig. 3. Dispersion curves taken at two different 
sample temperatures, The solid line corresponds 
to the curve taken at 2.03°K, while the broken 
line represents the dispersion curve observed at 
1.44°K, 


1.00 150 200 ar 


Fig. 3 shows the dispersion curve thus determined at 2.03°K. For comparison 
the curve observed at 1.44°K is included. The important feature of the new curve 
as regards the discussion of the data obtained at different helium temperatures 
is, that the minimum, x», occurs at approximately the same value at 2.03°K as 
was the case at 1.44°K. Within the experimental accuracy, and with due regard 
to the flatness of the curve in the neighbourhood of the minimum, one may regard 
the excitation energies determined from the experimental series at different tem- 
peratures as corresponding to the minimum. The error in the excitation energy, 
T°K, caused by this approximation is<0.1°K. 

Having establised this as a basis for discussion, we now proceed to the treat- 
ment of the data illustrated in Fig. 2. Table 2 contains the results of all the 
observations at the various helium temperatures. Equations (3) and (4) have been 
used to calculate line widths (dt), and 67’. 


Table 2. Line widths and excitation enrgies determined at an observation angle 
of 76.3° at different helium temperatures. 


: Flight : Ree 
oe time of | Observed Line aes aes Ax sch Roton Roton 
raat - scattered line width shift momentum| energy 
ok neutrons | (6¢)r usec | At psec (6 t)r OT ~ At UO 
ty sec psec SHA 
1.50 3305+ 6 | 100+ 7 236+ 8 1.89+0.04| 8.3+0.4 
1.65 3295+ 6 | 122+10 226+ 8 70+12 2.2+0.8 | 1.90+0.04] 8.0+0.4 
1.73 3289+ 6 | 142410 220+ 8 101412 3.2+0.6 |1.9020.04] 7.8+0.4 
1.85 3285+ 6 | 158410 | 216+ 8 122412 3.940.6 | 1.9040.04] 7.740.4 
1.94 3267+ 6] 200410 198+ 8 173412 5.7£0.5 |1.9040.04] 7.140.4 
1.98 325910 | 200415 190+ 12 173 + 16 5.70.7 |1.9140.05| 6.9+0.5 


2.03 | 3244410] 213+15 | 175+12 | 188+16 | 6.3+0.7]1.9140.05] 6440.5 
2.14 3197+10 | 298+15 | 128+12 | 281+16 | 9.8+0.7]1.92+0.05] 4.8+0.6 
2.17 3157+30 | 300+30 | 88+30 | 283+31 | 10.3+1.4 |1.9340.06] 3,341.3 
2.37 |(3140+15)| (330430) | (71416) | (314+31) | (11.6+1.3) (2.7+0.8) 
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Fig. 4. Variation of “gap’’ energy, A, as a function of the helium temperature. 


The roton energy values given in the last column of Table 2 thus very nearly 
correspond to the minimum energy A. This assumption is supported by the roton 
momenta belonging to the roton energy values as seen in the last-but-one column, 
which shows that the momentum values increase slightly with increasing helium 
temperature. But this is exactly what may be expected for the minimum accord- 
ing to the two complete measurements at 1.44 and 2.03°K. 


Discussion of the results 


According to Table 2 the minimum energy, A, necessary to create a roton 
decreases with increasing helium temperature first slowly up to about 1.8 or 
1.9°K, then more rapidly up to the j-point as shown in Fig. 4. Assuming that 
some unknown interaction between the excitations causes this change in A, it 
may be reasonable to suppose that the strength of this interaction is propor- 
tional to some power of the concentration, N,, of rotons per unit volume. Ac- 


os AR ep 
cording to the Landau theory, NV, VT-e *sT, Therefore we tried to fit our 
data with a function of the form A=4A- [i-a(vr-< r) i The result of such 


a calculation is that a very good fit to the experimental data is obtained with 
the following expressions: 


8.5 
is 85{1-75-1P-< "| for 1.4<7<1.9°K (5) 
Par — __85)186) 
8.5{1~ 1050. (Ve r) ge 18 <P<_2.15°K (6) 


A 
These expressions imply that the limiting value of = for very low temperatures 


B 
is 8.5°K, and that this value decreases linearly with N, as the temperature 
increases up to about 1.9°K. From this temperature and up to the A-point the 
temperature dependence is stronger and the correction term goes more like the 
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A point 


Sos 
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Fig. 5. Excitation line width as a function of helium temperature. The solid line corresponds to 
the fitted expressions (9) given in the text, while the broken line corresponds to the line width 
predicted by the Landau-Khalatnikov theory. 


square, N?, of the roton concentration. In the immediate neighbourhood of the 
A-point perhaps even higher powers of N, should be used. When the A-point is 
passed the scattered neutron distribution is represented by a wide gaussian dis- 
tribution (more correctly: the derivative of the scattered beryllium break is gauss- 
ian), the energy shift staying constant at a value of 2.7°K. 

The line width as given in the sixth column of Table 2 is represented graphi- 
cally in Fig. 5. According to the Landau-Khalatnikov theory, which assumes a 
constant A, the line width at least up to a helium temperature of 1.8°K may 
be calculated from 


h Pgoiiore 4.7 Wye 
cor eT, (7) 
Sun, 15 (22) n, A key 


oT 


where yw is the effective roton mass, 7, is the roton part of the viscosity, 9, is 


the roton part of the normal fluid density. Taking for 2 the value at 7’ =0°K 
B 


of 8.5°K, ~=0.16 mg. (see reference 2), © =1.90 A™* and 7,=1.1-10-° poise 
we find 
Yr yi 
6T~92-VT-e T (7') 


This calculated line width is shown as the lower curve of Fig. 5, the calculated 
values being about 4 of the experimentally observed ones. If one tries to fit 
8.5 


an expression of the type 67 =A: VT -e 7, using the asymptotic value of 8.5°K 
for A, one finds a good fit for A =302 up to a helium temperature of 2°K. From 
that temperature and up towards the A-point the line width increases more rap- 
idly. In accordance with the ideas discussed above, and assuming that some 
stronger interactions occur when the roton density increases to high values, we 


56 


ARKIV FOR FYSIK. Bd 15 nr 4 


—= | S0\a, 
tried to fit a function of the type A (vr =e r) to the observed line-width data. 
The final result of these fits, shown as the solid line through the experimental 
points in Fig. 5, is: 


8.5 
[ 302-VP-07 for T<2°K 


oT (8) 


8.5 


3 | 2060 - (vr Hr te for T>2°K 


There is thus a striking analogy between the variation of the line width and that 
of the minimum excitation energies with helium temperature. Up to the neigh- 
bourhood of 2°K both vary linearly with the roton density, N,, and for tempe- 
ratures above 1.8-2°K both vary more rapidly with temperature, indicating that 
some new phenomena come into play. Some sort of multiple interactions be- 
tween the neutron and two or more rotons simultaneously might perhaps not be 
excluded, a process analogous to multiphonon scattering of neutrons from crystals. 

The most striking result of the line-width determination is thus that the experi- 
mental values are about three times as large as the theoretically predicted ones. 
A consequence of this is, of course, that the mean free path of the rotons is 
three times shorter than the value calculated from the Landau-Khalatnikov the- 
ory. This discrepancy is too large to be explained by experimental errors, unless 
the method as a whole is masking large effects which should be interpreted dif- 
ferently. One such possibility would be that the roton line is asymmetric to- 
wards shorter wavelengths, resulting in too large apparent line widths and per- 
haps even too large line shifts. Within the given experimental resolution, however, 
no asymmetry of the line could be detected. Such a possibility can only be in- 
vestigated using a higher resolution and ideally by sending a narrow neutron 
“line” into the helium sample. 

Finally there is one point to be mentioned about the dispersion curve. An 
inspection of Fig. 3 shows that ‘there is a slight shift of the minimum towards 
higher momentum values simultaneously as the minimum energy A decreases for 
increasing helium temperature. A fit of a parabola to the excitation curve ob- 
served at 2.03°K gives for the Landau parameters A, and x: 


dash (6.7 +0.3) °K 
kp 
Mp = (1.94 +0.02) A~? 


fw =(0.13 £0.02) mye 


(9) 


A comparison of these results with the previous[2] data observed at 1.44°K 
A ‘ 

(2 =8.1°K, x,=1.90A", u=0.16 mi] shows, that the tendency for the disper- 
B 

sion curves when the helium temperature increases is, that the curves shift to- 

wards lower excitation energies, higher momentum values, and that the concavity 

of the parabolic region round the minimum increases. It should be pointed out 

that the previous [2] measurements at 1.44°K indeed showed that the parabola 

shape of the spectrum does not seem to hold except in a very narrow region, 

%) +0.2 A~*, round the minimum. This statement is confirmed by the fact that 
he data taken at 1.44°K are very well fitted by a function of the type 
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Fig. 6. Brueckner type dispersion curve fitted to the experimental data taken at a helium tem- 
perature of 1.44°K. 


r 
Toke. i Ve+2(2) sin 5% (10) 


(es Cy 


When c, = 4.74 +0.10, c, = 5.79 +0.06, cs = 2.047 +0.010. The result of such a fit 
is given in Fig. 6. The expression (10) is similar to the dispersion law given by 
Brueckner and Sawada except that a scaling factor included in c, has been ap- 
plied. It might be worth mentioning that this fitted curve predicts a maximum 
at 7’ =13.53°K and x=1.11 A“+, and that the slope of this curve at the origin 
is 19°K per A~', corresponding to a predicted velocity of first sound of 248 
metres/sec. This value thus surprisingly well agrees with the experimental value 
of 237 metres/sec. The full dispersion curve thus predicted agrees relatively well 
with the curve experimentally observed by the Los Alamos group [3], if the tem- 
perature variation is taken into account. 


APPENDIX 


A gaussian resolution function applied to a step function type spectrum 


Assume we have a gaussian resolution function of the type R(t—t?t’) and a 
neutron spectrum described by f(t). ¢ is the neutron time of flight and ?¢’ de- 


58 


ARKIV FOR FYSIK. Bd 15 nr 4 


f(t) 


R(t-t) 


' flight time 
% 


om --- 
° 


! 
\ 
Us 
t 


Fig. 1 (Appendix). Gaussian resolution function R(t—?’) applied to a triangular step function 
type spectrum. 


fines the centre line of the resolution function (Fig. 1). The intensity J (t’) ob- 
served at flight time t’ obviously is: 


I(t’) =[R(t—t') f(t) dt, (1) 
where R(t—t’) is normalized in such a way that 
fR@-t')dt=1 (2) 


_ (t=)? 
Ti we take for R(t—t’) the function c-e ®™ normalization gives: 


(t—t")® 


R(t-t')=—~=e ® 
9 a 
with 
a= th (4) 
2Vin2 


(Ot), is the full width at half maximum of the resolution function R. 
With these definitions we first calculate the shape of the observed spectrum 
for a step function neutron spectrum: 


(i) = ee for t<ty (5) 


1 for t>ty 


| The result according to eq. (1) for the observed intensity is then simply: 
| 


(t- ea 
Cla a (6) 
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Fig. 2 (Appendix). Spectrum J (¢’) resulting from the application of a gaussian resolution func- 

tion to the step function spectrum f(t). The half width of the resolution function, (0 ¢)), is very 

nearly given by the projection of the tangent at point of inflection on the ¢ axis. The point of 
inflection occurs exactly at to. 


This observed spectrum is illustrated in Fig. 2. As is easily proved, the projec- 
tion on the ¢ axis of the tangent at the point of inflection I(t’) defines a time 
interval which is 


= _ (bt) Le : 
aVa= 5 ind 1-08 (580 (7) 
For this case we thus find that the point of inflection of J (t’) defines the “‘cut- 
off’ t,; exactly and the projection of the tangent on the ¢ axis very nearly defines 
the full width (d#)) at half maximum of the resolution function. 

Next we suppose the spectrum f(t) is triangular as in Fig. | and is given by: 


t—t 
1 for th<t<t, 


f=, 4-b (8) 


0 for t<t, and t>t, 


The integration according to eq. (1) then results in intensity distributions I (t’) 
of the shape illustrated in Fig. 3. These curves are characterized by two facts: 
a) The projection on the ¢ axis of the tangent at the point of inflection still 
very well defines the resolution width (dt), and particularly so when (6 t))><t, — ty. 
b) The inflection point no longer corresponds exactly to t, but is shifted to- 
wards a smaller ¢ value, t as shown in Fig. 3. This error in the definition of 
t), denoted by At in Fig. 3, increases as t,—t, decreases for a constant (dt)y but 
remains small as long as (dt))><t,—t. Fig. 4 gives the error At as a function of 
ae As will be seen from this figure the correction in the actual case, when 
0 
the spectrum f(t) is well approximated by eq. (8), with (df))~+100 usec and 
t, —t) ~ 1000 usec, is only 3—4 usec. 
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Fig. 3 (Appendix). Spectrum J (¢’) resulting from the application of a gaussian resolution func- 
tion to the triangular spectrum /(t). The observed “break” position, f), is now shifted by an 


amount At from the true “break” position. 
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ig. 4 (Appendix). The correction, At, to be added to the observed “‘break’’ position as a func- 
t,—ty 


(Ot) 


tion of the ratio 
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As a next step in our analysis we let our spectrum f(t) have a gaussian 
shaped front: 


1 EEO 
(== fe on dt” 9 
bVa e 

to 
The quantity b is according to the preceding discussion, eq. (4), defined by 
(dt); 

a pr 4’ 
2Vin2 ee 


where (dt), is the full width at half maximum of the gaussian line shape now 
assumed to have distorted the step function spectrum f(t) to the shape illustrated 
in Fig. 2. If we now observe the spectrum f(é) given by eq. (9) with a spectro- 
meter having a resolution function R(f—t’) as given by eq. (3) we obviously find 
for the observed intensity: 

(t'—)* 


Lh ee eh atfePae (10 


—0 t,—t 


where «=t’’—t. As is easily shown, this integral may be transformed so that the 
observed intensity results from: 
1 BS a 
LG) eee fe wrerat 11 
Vx (a? + 8) ( 


But this means that the full width at half maximum of the observed gaussian 


line (which is at ' is proportional to Va?+b. By using the definitions of a 


and 6 we find 
(8 t)ovs = V(b t)o + (6t)i, (12 


which is just the well known addition theorem for gaussian functions. For such 
a spectrum as f(t) given by eq. (9) and reproduced by our apparatus as given 
by eq. (11) we may still use the procedure outlined in the preceding section tc 
find (dt),ps, and from that we may calculate the line shape defined by (61),, if 
(dt), is known. Also the point of inflection of I(t’) defines f, exactly. 

When the spectrum f(t) is not given by eq. (9) but has a triangular shape 
rounded off by a gaussian shaped front the same conclusions remain valid, ex- 
cept that the point of inflection now gives a value t) which should be corrected 
as given in Fig. 4 to give the true “cut-off” position ft). 


ACKNOWLEDGEMENTS 


The authors are very much indebted to Dr. Harry Palevsky, Brookhaven, U.S.A, for stimu: 
lating discussions. We have also profited from a fruitful correspondence with Professors Landau 
and Khalatnikov, USSR. 


AB Atomenergi, Department of Physics, Stockholm. 


62 


ARKIV FOR FYSIK. Bd 15 nr 4 
REFERENCES 


[1] H. Patevsky, K. Ornes, K. E. Larsson, R. Pautt and R. StepMan: Phys. Rev. 108, 
1346 (1957). 


[2] H. Patevsxy, K. Ornss, and K. E. Larsson: Phys. Rev. 112, 11 (1958). 

[3] J. L. Yarnecn, G. P. ARNoxp, P. J. Benpt and E. C. Kerr: ‘“‘Energy vs. Momentum Rela- 
tion for the Excitations in Liquid Helium”. Report at Kammerlingh Onnes Conference 
on Low Temperature Physics in Leiden, June 23rd—28th, 1958. 

[4] M. Conen and R. P. Feynman: Phys. Rev. 107, 13 (1957). 

[5] L. Lanpav: J. Phys. USSR, 5, 71 (1941). 

[6] R. P. Feynman: Phys. Rev. 94, 262 (1954). 

|7] K. A. BRuECKNER and K. Sawapa: Phys. Rev. 106, 1117 and 1128 (1957). 

[8] K. E. Larsson, R. Stepman and H. Patrysxy: J. Nuclear Energy, 6, 222 (1958). 


Tryckt den 27 januari 1959 


Uppsala 1959. Almqvist & Wiksells Boktryckeri AB 


63 


erty ica 
AS) SUE 7 
:( hd abe eae oe 


: z JER Baggs ake 


